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Abstract 



In the exact propagator for finite temperature gluons the location of the trans- 
verse and longitudinal poles in the gluon propagator are unknown functions 
of wave vector: wt(&) and coL(k). The residues of the poles, also unknown, 
fix the normalization of the one gluon vector potential and thus of the field 
strength. The naive energy density E ■ D + B ■ H is not correct because 
of dispersion. By keeping the modulations due to the source currents the 
energy density is shown to be ojt/V and lol/V regardless of the functional 
form of u>T{k) and u)T,{k). The momentum density is k/V . The resulting 
energy-momentum tensor is not symmetric. 



I. INTRODUCTION 

Zero temperature. At zero temperature the canonical operator formalism guarantees that 
a gluon of frequency u has a field configuration with energy density u/V. The same result 
can be obtained from the propagator because the residue of the pole in the propagator at 
uj = k fixes the potential of a transverse gluon to be 

A = e(^) l/2 e ^r-kt) (L1) 
V2JfeW v ; 

The resulting electric and magnetic fields have an energy density {\E\ 2 + \B\ 2 )/Z 3 = ui/V. 
In a covariant gauge the factor I/Z3 may be viewed as the constant value of the electric 
permittivity and magnetic permeability 1/e = \i = Z 3 for on-shell gluons in the vacuum, 
and the energy density is e\E\ 2 + \B\ 2 /fi = uj/V. It is important to realize that Z3 is a 
complicated function of the coupling constant and depends on the infrared and ultraviolet 
regularization. But the statement that the energy density is u/V is exact. It is useful to 
restate this result. Suppose the full zero-temperature effective action T(A,ip) were known 
and from this the energy-momentum tensor T^^A^ip) could be constructed. The energy 
density T 00 would contain terms quadratic in A, as well as cubic and higher powers and also 
terms containing ip. If the value of A is put at the one-gluon value (1.1) it is the quadratic 
term in the exact T 00 that would equal u/V. Terms in T 00 that contain higher powers of A 
do not describe the energy of a single on-shell gluon and should not be included. 

Hard thermal loop approximation. At non-zero temperature the analysis of the field en- 
ergy of thermal gluon oscillations has previously been done within the hard thermal loop 



approximation of Braaten and Pisarski [1-3]. This is a consistent high-temperature ap- 
proximation to QCD in which only the order T 2 contributions of all one-loop diagrams are 
retained. The gluon self-energy is gauge invariant and satisfies K^W^ = 0. The poles in 
the gluon propagator are not on the light cone. For the two polarizations transverse to 
the wave- vector k the pole is at a frequency uri^k) that is proportional to gT but with a 
complicated dependence on k. There is another pole in the propagator for gluons polarized 
longitudinally, i.e. parallel to k, and this pole is at a different frequency UL(k). The energy 
density of the fields is not proportional to \E\ 2 + \B\ 2 . The analysis of the energy and mo- 
mentum has been specific to the detailed properties of the hard thermal loop approximation. 
Because the off-shell amplitudes with any number of external gluons have the remarkable 
property of being gauge-fixing independent and obeying abelian Ward identities [2,4], it 
is possible to construct effective actions that generate all the hard thermal loops [3,5-11]. 
The effective actions are highly non-local and consequently the canonical procedures for 
constructing the energy-momentum tensor do not apply. The full energy- momentum tensor 
contains quadratic and all higher all powers of the vector potential since they summarize 
all higher correlations. From the Braaten-Pisarski form of the action [3], I constructed a 
conserved energy momentum tensor [12] but the quadratic part of that T 00 does not give the 
physical energy density u a /V for the two modes. Other authors have constructed different 
conserved energy-momentum tensors [13-15], also containing all higher powers of A^, and 
have generalized them to curved space-time [16]. The non-uniqueness of T^ u is expected 
since a term d a X aixv will be conserved with respect to <9 M for any X afil/ that is antisymmet- 
ric under interchange of a and fi. Blaizot and Iancu [15] found an energy- momentum tensor 
that does give the correct energy densities for the transverse and longitudinal plane waves. 
It differs from [12] by such a total derivative term. In these analyses the detailed struc- 
ture of the quadratic, cubic, and all higher correlations are related by the Ward identities. 
Consequently even in the quadratic approximation the specific k- dependence of uxik) and 
ujL(k) is essential. 

The purpose of the present paper is to go beyond the hard thermal loop approximation in 
computing the field energy of the thermal gluon oscillations. This question is simple enough 
that it can be answered exactly just from the structure of the exact gluon propagator [17]. 
The inverse of the propagator gives the quadratic term in effective action. As emphasized 
above, experience with the hard thermal loop approximation shows that the spatial and 
temporal non-localities of the quadratic action prevent a unique construction of even the 
quadratic terms in energy-momentum tensor. In this paper a physical argument will be used 
along the following lines. First, the quadratic terms in the action are used to obtain the 
linear field equation 5T(A)/5A IM (x) = J^(x) for the gluon potentials in the presence of an 
arbitrary color current J M . In coordinate space this will turn out to be Maxwell's equation 

d,H^ = J£ ee + 0{J 2 ) (1.2) 

in a media characterized by an electric permittivity t{K) and a magnetic permeability fi(K) 
that are related to the gluon self-energy. Thus even in the exact theory, the gluon oscillations 
have many features in common with the classical electrodynamics of a classical plasma. 

— * — * 

Although the color electric field E and the color magnetic induction field B are produced 
by all the quarks and gluons, the contribution of the moving quarks and gluons to the self- 
energy tensor produces polarization and magnetization currents which are incorporated into 
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the electric displacement field Di = e(K)Ei(K) = Hqi and the magnetic field He(K) = 
Bi(K)/ fi(K) = e£ mn H mn /2. The sources of these macroscopic fields is referred to as the free 
current J f ^ ee . The label 'free' does not mean trivial, but rather that the polarization and 
magnetization currents are not included in this current as that would be double counting. 
In terms of Feynman diagrams J£ ee is one-gluon-irreducible because the self-energy effects 
of the gluon are contained in the electric permittivity e(K) and the magnetic permeability 

/'(/<)• 

The next step is to construct the energy and momentum density of the gluon fields. The 
energy and momentum of a source-free plane wave cannot be simply identified. Instead 
it is necesary to examine field configurations that are not completely source-free and thus 
are not precisely plane waves. The energy and momentum of such modulated waves can 
be determined from the currents that radiate those fields. The essential physical step is to 
realize that the power supplied to the fields and the force exerted on the fields is given by 
the real part of J^l e F a v for harmonic time dependence. This is expressed in the identity 

d ^ H ; x F^) + \mjerF^) = r airee F™. (1.3) 

Since the right hand side gives the rate of energy and momentum input to the fields by the 
currents, the left hand side must describe the energy and momentum content of the field 
configurations. To use (1.3) one must solve Maxwell's equations in the presence of a weak 
source J f " ee that supplies energy and momentum to the fields. The source produces modu- 
lated fields of the form exp(ikz — iut) f(z, t) where E, B, D, H all have different modulation 
functions / because of the dispersion. Keeping terms to first order in df/dt and df/dz 
reduces the left side of (1.3) to a total derivative: 

^ = Re(J: free F-) (1.4) 

Since there is absorption in the plasma the energy and momentum of any traveling wave will 
not truly be constant in time. In applying (1.4) it is necessary to neglect the damping or, 
equivalently, to neglect the imaginary parts of e(K) and n(K). This is a physical limitation. 
Thus although the poles in the gluon propagator occur at complex values of frequency, only 
the real parts uixik) and uJiik) will be retained. In the limit that the source J f " ee — > 
the modulated fields become plane wave fields whose magnitude is completely fixed by the 
residues of the poles in the propagator. For these plane wave modes of single gluons the 
energy-momentum tensor in either mode is 

a = T or L (1.5) 

where vt = dux/dk and vl = du L /dk are the group velocities. By the theorem of Kobes, 
Kunstatter and Rebhan [18] the functions u>T{k) and ujx(k) are gauge invariant. The com- 
ponents of (1.4) satisfy simple physical checks. First, the energy density T 00 = uj/V and the 
momentum density T 0n = k n /V. Because the energy and momentum are transported at the 
group velocity, T m0 = v m T 00 and T mn = v m T 0n . No assumption is made about the form of 
the dispersion relations ur(k) and UL(k). (Note that T^ v would be symmetric for dispersion 
relations of the form form u = (k 2 + m 2 ) 1//2 because then vuo — k. That is certainly not the 
case here.) 
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The analysis performed here is intended to be simple and physical. It is a generalization of 
old results in optics and classical plasmas [19-21]. In those contexts e is assumed independent 
of k and /i — 1. Most treatments are rather unphysical in that they omit discussion of the 
external current that supplies the energy and momentum content of the fields even thought 
the modulated fields are only possible with external sources. 

The remainder of the paper is organized so that all physics is deduced from the exact 
thermal gluon propagator. Sec 2 uses the poles in the propagator to obtain the normalized 
transverse and longitudinal plane waves. Sec 3 derives Maxwell's equations from the gluon 
propagator and the relation of e(K) and /jl(K) to the gluon self-energy. Sec 4 deduces the 
energy and momentum for the transverse mode. Sec 5 does the same for the longitudinal 
mode. 



II. POLES IN THE THERMAL GLUON PROPAGATOR 

A. Gluon propagator in general linear gauges 

One of the first obstacles in constructing the exact propagator for finite temperature 
gluons is that K f J\^ v ^ [22-25]. It was shown in [17] that for any linear gauge- fixing 
condition the Slavnov- Taylor implies that Det[ — K 2 b il v + K^K V + IPJ = 0. Local gauge 
invariance requires one of the four eigenvalues to be zero for all . This provides a non- 
linear condition on the self-energy tensor [24-28]. It is easiest to work with a gauge- fixing 
that preserves the rotational invariance of the plasma rest frame. These include covariant 
gauges, Coulomb gauges, and temporal gauges. In rotationally covariant gauges the gluon 
polarizations that are transverse to the spatial momentum k are isolated by the projection 
tensor 



/0 \ 

rest ~~ I -5 mn + k m k n J ' ' ' ' 



The gluon four-momentum is and a related vector in which the energy and momentum 
components are reversed is denoted by K^: 

K» = (u,k) K» = (k,uk) (2.2) 

The new vector satisfies K 2 = —K 2 and K • K — 0. The same construction applies to any 
vector whose spatial component is along k: 

V» = (v ,vk) V» = (v,v k) (2.3) 

These satisfy V 2 = —V 2 and V ■ V — 0. For a linear gauge-fixing condition that preserves 
the rotational invariance the inverse of the free propagator has the form 

BtaW] =-K 2 ! r + K»ir-— (2.4) 

where the gauge-fixing vector has spatial components parallel to k. The free propagator 
is 
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j^xv pup" K^K V 
DZe(K) = + (F . K)2 - t {F . K y ( 2 - 5 ) 

Familiar examples of this include covariant gauges [F^ = K^}, Coulomb gauges [F M = 
(k, 0,0,0)], and temporal gauges [F M = (0, k)]. Including interactions gives a proper self- 
energy tensor of the form 

w v = u t a^ + iu-^- + n c — + n D -^- (2.6) 

because of rotational invariance [24-28]. (For QED and also for QCD in the hard ther- 
mal loop approximation He = Hd = 0.) The Slavnov- Taylor condition leads to He = 



a J (K 2 — Ul)TId- and the full propagator is 



where = K^y/K 2 — Hi + oK^\JH d . See also [22,27]. The second denominator may be 
written directly in terms of the self-energy tensor as 



(F-H) 



2 



(F ■ Kf + F^IP" (2.8) 



K 2 

The longitudinal denominator for three simple gauge choices are as follows: 

Covariant : ^ = (K 2 ) 2 + K^K U W U (2.9) 

Coulomb : (F "f )2 = k 2 {k 2 + IT 00 ) (2.10) 

Temporal : ^ = k o + k m k n TT n (2.11) 

It is sometimes convenient to parametrize the gauge-fixing vector as F^ = fi K^ 1 + / 2 so 
that F ■ H = K 2 (fWK 2 - n L - af 2y /TlE). 



B. Transverse plane waves 

The poles in the propagator ( |2.7| ) determine the freely propagating plane wave solutions. 
For the transverse mode the value of u for which K 2 = Ut defines the complex energy 
uj t = ujt — ijT- At this pole 

w-u4: D mn (K)^ . Zt - (S mn -k m k n ), (2.12) 

2uj t (uj — lu t ) 



where the wave function renormalization factor is defined by 

(2.13) 



2"* _n a 
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The residue of the pole fixes the normalization of the vector potential for one transverse 
gluon to be 

A = -ie(-%) 1/2 e^-^ (2.14) 

with V the quantization volume. The phase factor % has been chosen for later convenience. 
As usual e is a unit vector satisfying k • e — 0. For definiteness the direction of propagation 
will be chosen as the z axis and the direction of polarization as the x axis: 

Ax= dJ^\ xt2 ^-^) (2.15) 



The corresponding fields are 



Ex = u cjJ^) 1/2 e <V"-**Q By = k{^-) 1 ' 2 e «<*-«4*> (2.16) 

It will be shown later that these fields propagate freely in that they satisfy Maxwell's equa- 
tions with an electric permittivity e and a magnetic permeability \x but with with no free 
current. 



C. Longitudinal plane waves 



When [F ■ H) 2 = the propagator (2/7 ) has a different pole which will be called lon- 
gitudinal because the resulting electric field is along the spatial wave vector k and there is 
no magnetic field. The value of to at which the denominator (F ■ H) 2 vanishes defines the 
complex energy oj c L) which has real and imaginary parts uj c l = ujl — ijL. The propagator at 
this pole behaves as 



CO 



D^(K) 



2u c l (uj - u c L ) 



K 2 



(2.17) 



where the £ dependent tensor in (|2.7f ) is omitted. The wave function renormalization factor 
for the longitudinal mode is 



2uj c , 



d /F ■ H\2 

aJv k 2 



(2.18) 



The coefficient of the pole fixes the normalization of the vector potential for one longitudinal 
gluon to be 



1/2 



pV e -i(k-r-uj c L t) 



\2ulVK 2 Lj 

where the phase factor — i has been chosen for later convenience. Since 



(2.19) 



hK» + f 2 K^ 



(2.20) 
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the potential is 



f (7 K 2 \ 1 ^ 2 

A* = i J M^£) e~ iKL - x + (•••) K». (2.21) 

K \ 2jLO V J 

The omitted terms proportional to do not contribute to the abelian field strength. Only 
the potential along w M is non-trivial and in the plasma rest frame the only non-vanishing 
field is 

^ = /i(f^) 1/! 1 j {kz - lt) - (2-22) 

Since there is no magnetic field, the way in which this mode transports momentum is par- 
ticularly subtle. 



III. MAXWELL'S EQUATIONS FROM THE GLUON PROPAGATOR 

The previous section displayed the plane waves that come from the propagator poles. 
This section will show that the off-shell gluon propagator implies Maxwell's equations in a 
media. 



A. Electric and magnetic phenomena 

Current conservation. The gluon propagator determines the quadratic term in the gen- 
erating functional for one-particle-irreducible diagrams in the presence of a colored source 
current J": 

W(J) = \J d'xd'y j;{x)D^{x - y)J a v {y) + 0(J 3 ). (3.1) 

All four components of J M are independent variables, not constrained by current conserva- 
tion. In order to obtain the equation of motion for the gauge potentials it is necessary to 
impose color conservation, D a ^J c ^ = 0, where 

D°f = 5 ac d„ - gf abc - I \ . (3.2) 

To impose the constraint it is convenient to introduce a Lagrange multiplier function if)(x) 
and vary all four components of J M independently in the functional 

W'(J) = W(J) ~\jd 4 x i,{x){D»J,) a {D»J^ (3.3) 
The resulting vector potential = 5W /SJ^ is 

A a »(x) = J d 4 y D^(x - y)J a v {y) + ^{x)d^d v J*{x) + 0{J 2 ) (3.4) 
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In momentum space 



A a »{K) = D^(K) - ^K»K V J a v {K) + 0{J 2 



(3.5) 



where \I/ = ip(—K). From now on the color indices will be omitted. To solve for the Lagrange 
multiplier function, contract both sides of the above with the vector and use the fact 
that Q satisfies F^ v = -£H»J U /(F ■ H) to get: 



H"J V (K) -*(F-K) K V J U (K) 



F-H 

Substituting this into (|3~5| ) and rearranging gives 



A'"(K) = D^(K)+£ 



(F.K)(F-H) 



MK), 



where 



A ,fi (K) = 



K»F V 
K-F 



A V {K). 



(3.6) 



(3.7) 



(3.8) 



This potential satisfies the gauge condition F^A'^^K) = 0. In coordinate space ( |3.7| ) is the 
integral equation for A' 11 . By multiplying from the left by the tensor shown below, ( p.7|) can 
be converted into the differential equation 



(-K 2 + n T )A* + ( T =-) ——}A' fl (K) = JgjK) 



f t K 2 ) K 2 



(3.9) 



where the free current is defined as 



J^(K) = J a {K)+£Jl D 



k a k v j v [K) 
~TKk 



2\2 



(3.10) 



and satisfies K a J^. ee (K) = 0. The transverse plane waves of Sec 2.2 at K 2 — Ut = 
propagate with Jf ree = 0. Similarly the longitudinal plane waves of Sec 2.3 at (F ■ H) 2 = 
propagate with Jf ree = 0. It is sometimes useful to have the inverse form of (|3.9|) : 



A'(K) 



A 



fJ,U 



K 2 -Tl T 



+ K 2 



F F 



(F ■ Hf 



(3.11) 



Identification of e and /j,. Equations ( |3.9|) are the inhomogeneous Maxwell equations in 
a media characterized by an electric permittivity e and a magnetic permeability fi. To see 
this one needs to separate the electric and magnetic component of K^A' (K) by using 



k K^A'^(K)(K) = iE\K) + u{8 lm - k e k rn )A' rn {K), 
where E(K) = -ikA'°(K) + iul'(K)(K). The a = component of O is 



(3.12) 



e(k,u) ik.E(K) = J° ee , 



(3.13) 
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with a gauge- dependent electric permitivity given by 

1 / F ■ H\ 2 F F 17^ 

The electric displacement field is D(K) = eE(K). Next comes the Ampere-Maxwell equa- 
tion, which is contained in the spatial components of ( |3.9j ): 

(K 2 - Il T )(-S £m + k r k m )A' m {K) + e(k,Lu)ujk e K»A'^K) = j/ ree . (3.15) 

Using (|3.12| ) this becomes 

1 



fi{k,u] 



ik x B(K) + e(k, uj)iuE(K) = J cxt , (3.16) 



where B(K) = ikxA'(K) is the magnetic induction field and the gauge-dependent magnetic 
permeability is 

1 1 [-K 2 + U T + u 2 e(k,u)}. (3.17) 



fj,(k,u) k 2 
The magnetic field is H(K) = B(K)//i 



B. Energy and Momentum Conservation 

In coordinate space the inhomogeneous equations ( |3.13 ) and (|3.16| ) take the familiar form 

■/free- (3.18) 



V • D = J° ee 



The two homogeneous equations, V • B — and V x E + dB/dt = 0, are true for any A^. 
These results can be expressed in terms of the tensors 



( 


—D x 


~Dy 


-D z \ 


D x 





-H z 


Hy 


Dy 


H z 









~Hy 


H x 


o) 



I 


—E x 


~Ey 


-E z \ 


E x 





-B z 


By 


Ey 


B z 





— B x 


\E Z 


-By 


B x 


0/ 



(3.19) 



Maxwell's equations in the rest frame of the plasma can be written 

ga pPl _|_ QP _|_ g-y paf3 _ 



Tr 

° free 



(3.20) 



In order to identify the energy and momentum in the fields, the essential physical step is 
to recognize that the work done on the fields by the currents and the force exerted on the 
fields by the currents are given by the real part of J* { ree F ai/ for harmonic time dependence. 
This is expressed in the identity 



d ^ H ; x F x n 



l -H{ a {d»F^) 



J* fp au 



(3.21) 
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which follows directly from ( |3.20|) . The v = component is 



V • {H* x E) 



dD* 
~dt 



■E 



H ■ IK 



free ' E- 



(3.22) 



The right hand side is the power supplied to the fields by the color charges in the source 
current. The left hand side contains the spatial and temporal changes in the field energy as 
will later be shown. The v 



d 



3i {D " X B 



1, 2, 3 components of ( 3.21 



_ T°* F — 



arc 



D*VE j 



J free X 



(3.23) 



The right hand side is the force acting on the fields by the color charges. This external force 
is responsible for changes in the field momentum that are contained on the left hand side. 



IV. TRANSVERSE GLUON OSCILLATIONS 
A. Modulated fields radiated from a current 



The transverse plane waves ( 2.16|) that correspond to poles in the propagator at K 2 = Ut 



are exact solutions of Maxwell's equations ( |3.9D with J f ^ ec = 0. As discussed in Sec 1, to 
identify the energy and momentum it is necessary to consider a modulated wave of the form 

Ex = e i(kz- ur t) E T {z,t), (4.1) 

where the amplitude function ET(z,t) is almost constant. The only significant variation in 
the amplitude function should be when the wave is emitted or absorbed (i.e. at large \z\ 
and large \t\). In Fourier space the amplitude is 

E T {z,t) = J dad(3 C T {a,f3) e i{az ~ pt \ (4.2) 

where Cr(a,/3) should have support only for very small wave vector a and very small 
frequency (3. Note that this is not a superposition of source-free fields, which would be the 
case if k + a produced a pole at ujt + P- Here a and (3 are independent variables; the wave is 
slightly off the mass-shell. The total wave vector and total frequency in ( f4.1|) will be denoted 

k' = k + a to' = cot + P- (4.3) 

Faraday's law requires that the magnetic induction field be 

B v = e i{kz ~ UTt) ( dad{3 —C T (a,p) e i{az ~ pt) . (4.4) 
J U) 

The constituitive relations then fix the electric displacement field and the magnetic field to 
be 

D x = e W*-<>**) f dad(3 e{k' , u')C T {a, p) e i{az ~ pt) (4.5) 
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Hy = < 



" / — / dad(3 ^ ?C T (a,/3) e^-M. 

uj'fi(k', uj') 



(4.6) 



The current which radiates these modulated fields can be easily computed from (|3.18| ). The 
only non- vanishing component of J f ^ ee is J x with 



J, 



x free 



,12 



Akz^ T t) dadf3 -u'e(k',uj') C T (a,p) e i{az ~ pt) . 



(4.7) 



In the plane wave limit CV(a,/3) — > 5(a)S(j3)E and the source will vanish, i.e. J x f re e — > 0. 

The fields E x , D x , B y , H y each have slightly different modulations multiplying the 
exp(ikz — iujTt) factor. It is these slight differences that will account for the energy balance. 
For example, in the integrand of B y the factor k' /uj' can be expanded to first order as 



By = r 



M/ '- ' 1 / dad(3 (— + — a - \(3 + • • -)C T (a, (3) e i( - az - 0t) . 

Ijjrp Ijjrp UJrp 



(4.8) 



Only these terms are important since Cr(a,/3) has support only for very small values of a 

(4.9) 



and (3. Using flOf ) this can be written as 

By 

In the same fashion 



i(kz-u T t)(]^ ET _ J^dEr _ %k_dEr_ + 
^ujt ujt dz uo\ dt 



« kz -^(eE T -i^^ + i^^ + . 



Hy = i 



(- 



k 



-E T — i 



dk dz du dt 



d r k \ dE T , d / k \ dE T 
Ok^u/i' dz dkAu;//' dt 



+ 



(4.10) 



(4.11) 



Note that e, /x, and the partial derivatives are to be evaluated at to = oj t . If e and \i were 
constants, then the above relations would reduce to D x = eE x and H y = B y j ji. 



B. Energy Conservation 

For the TEM mode with E polarized along the x axis, the work-energy theorem (|3.22|) 



is 

d dil* dR 

Y MH;E X ) + M^fE x ) + MH* y ^f) = -Re(J* xfiee E x ) (4.12) 

As discussed in Sec 1, it is necessary to treat e and /i as real in order to ultimately obtain 
a time-independent energy. For harmonic fields and currents the real part of each product 
gives one-half the instantaneous value as will be seen below. To first order in the derivatives 
of Et, the spatial divergence term gives 

sS'W = ir T h E ^> (413) 
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as expected. However the time derivatives in (|4.12|) become 



dD* 



d 



d \E T \ 2 



d^ [UJe) W 



— UJ 



dt d \E T \ 
dk~dz~^2 



(4.14) 



_d_(ti\ d \E T \ 2 



k 2 f d 1\ d \E T \ 2 
+ — ( — -)■ 



duj^jiuj' dt 2 uo T ^dk \i> dz 2 



(4.15) 



It is because of spatial dispersion, specifically because e and /i depend on k, that the right 
hand sides of Q4.14D and (|4.15| ) contain spatial derivatives of \E T \ 2 . The sum of ( |4.13|) - (f4.15|) 



is 



2/.- ^de_ 



k 2 d li d \E T \ 2 



Imjt dk ljt dk fi^ dz 2 



+ 



d , s d f k 2 ^d \E T 12 
— (wei ' — 

OUJ 



du^fjw^dt 2 



J T\ 



It is convenient to define f(k,u) = cue — k 2 /(ufi). Then the work-energy theorem can be 
stated as 



df d \E T \ 2 df d \E T \ 2 
~ dk dz 2 + du dt 2 

The energy flux and energy density are 



— Re (.£/£</* 



T 30 (z,t) = - 



dk 



Ef\ 2 



T™(z,t) = ^-\E T \ 2 . 



(4.16) 



(4.17) 



Relation (|3.17| ) allows / to be expressed entirely in terms of the transverse self-energy: 
f(k,u) = [K 2 — ELr)/a>. The k and u derivatives in ( |4.16[ ) are to be evaluated at ut, i.e. at 
/ = 0. Therefore the energy is transported at the group velocity of the wave: 



dixlT 
~~dk 



(4.18) 



Because df / duo is related to the residue Z T through ( p,13| ) , the energy density can be written 

(4.19) 



T 00 (z,t) = ^\E T \ 2 . 



2 



Plane Wave Limit. The final test of T 00 is that it give the correct energy density in 
the plane wave limit: J^ t — * 0. The electric field magnitude of the plane wave TEM gluon 
oscillation is fixed by the propagator to be \Et\ 2 = ujtZt/2V in ( |2.16| ). Thus the energy 
density of one transverse gluon is T 00 = ujt/V as it should be. 

C. Momentum Conservation 



For the TEM mode the equation for momentum conservation ( 3.23 ) is 

,dE x 



d (n*R) + dH *yn 

Wt {D x B y ) + —D y 



D* 



x dz 



T* R 



(4.20) 
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Substituting the fields ( |4.9|) - ([4.11|) and taking the real part gives 



k df d \E T \ 2 k df d\E T \ 2 _ 
~ ~ T Tk d~ Z — + ~ T d^ ~ -MJ**~Pv)- (4-21) 

The momentum flux and momentum density are therefore 

T 33 = - — %\E T \ 2 T m = —^\E T \ 2 (4.22) 

LUt Ok LOT OLO 

The momentum is transported at the group velocity of the wave: 

r3 W^) = -(|/IU = ^. (4.23) 



Using (|2.13| ) one can write the momentum density as 



T 03 = Jk_ | St |2 (4>24) 

With the plane wave value of Et this gives the correct value T 03 = k/V . 

V. LONGITUDINAL GLUONIC OSCILLATIONS 

A. Modulated fields radiated from a current 

The longitudinal plane wave (2.22) comes from a pole in the propagator at (F ■ H) 2 = 
satisfying Maxwell's equations ( |3.9| ) with Jf ree = 0. To identify the energy and momentum 
content one needs a modulated wave 

Ez = e i{*z-"it) E L (z,t), (5.1) 

where the amplitude E^z^) is slowly varying in space and in time. The only significant 
variation should be at \z\ — > oo and \t\ — > oo. It has the Fourier representation 

E L (z,t) = J dad(3 C L {a,[3) e t{az - pt \ (5.2) 

where the function Cl{ol, /3) has support only for small a and small (3. The electric field has 
an associated displacement field 

D z = e l{kz ~ ULt) [ dadp e(k', uj')C L {a, (3) e l{az ~ pt \ (5.3) 



where k' = k + a and cu' = cul + (3 as before. Both H and B are identically zero. The 
modulated wave requires an external gluon current with components 

Jzbee = e l{kz ^ Lt) f dad(3 iu'e(k', to')C L {a, (3) e^ az ^ (5.4) 
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Jo free = e i{kz ^ Li) f dad(3 ik'e{k', Lu')C L (a, (3) e l{az - pt) (5.5) 



In the plane wave limit Cl — > S(a)S(j3) and therefore e(k',u>') — >■ e{k,u)L) = 0. This means 
that 7J 2 — ► 0, J 2 — * 0, and Jo — > eventually. 

Before taking the plane wave limit, the modulations in D z are different than those of 
E z because the former includes the field of the polarization charges. To approximate D z , 
expand e(k',u') in a Taylor series: 

D z = e -(**-<--.-) [ dad(5(a^- + (3^- + • • >\C L (a,l3) e i{az - pt) . (5.6) 
J v ok dco ' 

In terms of El this means that 

The derivatives of e are evaluated at e = 0. 



B. Energy Conservation 



Since the only nonvanishing fields for the modulated longitudinal mode are E z and D z , 
the work-energy theorem (|3.22|) is 



-- ,dD* 
Re{—^E Z 
at 



-Re(J* free E z ). 



(5.8) 



Despite its rather trivial appearance this equation does describe the flow of energy in the 
plasma. Because of 



dt 



de dE* L de dE* L 



' dk dz du dt 
Taking the real part, with e assumed real, gives 



de d \E L \ 



de d \E L \ 



'dkdz 2 "dudt 2 
Thus the energy flux and the energy density are 



-Re(E z J* hee ). 



T 



30 



96 I TP I 



de 

1 — ^LTT- \El\ 

duj 



(5.9) 



(5.10) 



(5.11) 



with the derivatives evaluated at e = 0. The energy is transported at the group velocity of 
the longitudinal wave: 



T™(z,t) T™(z,t) = - ^ 



de I de \ cLol 
dk' du^=o ~ ~dk 



Because of (2.18) and <KW) 



de_ 

duj 



e=0 



1 2u L 
fiK 2 Z L ' 
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(5.12) 



(5.13) 



The energy density is therefore 



2lu 



2 



T"' = jfgfc I^P (5.M) 

Plane Wave Limit. Whenever the external current is present the space and time deriva- 
tives of \EL(z,t)\ 2 prescribe the energy transport within the plasma. When the external 
current is removed, the value of T 00 remains constant. The value of that constant is the 
energy contained in the wave. The pole in the propagator fixes the amplitude for the free 
longitudinal gluon mode to be \E L \ 2 = f 2 Z L K\j \2u L V) as shown in (|2.22| ). The energy 
density of this mode is therefore T 00 = uj^/V. 

C. Momentum conservation 

For the longitudinal mode the force equation (|3.23|) is 

f)D* 

-df E * = - J ^E z . (5.15) 

Note that although J z f ree is non-zero, it exerts no force on the plasma charges in this mode 
because B = 0. Substituting (|5.7|) and taking the real part gives 



t* + »l^!! = _Re( j£>ft ). (5.16) 



dk dz 2 du dt 2 
Thus the momentum flux and momentum density are 



T 33 = -k^- \E L \ 2 T 0S = k^\E L \ 2 . (5.17) 
Momentum is transported at the group velocity of the wave 

^/^--(l/*)^. (5,8, 
The momentum density can also be written 

In the plane wave limit \E L \ 2 = flZ L K 2 L /{2u L V) so that T 03 = k/V, which is the correct 
result. 



VI. CONCLUSION 

Regardless of the functional dependence of u>t and u>l on momentum, temperature, and 
quark masses the energy density carried by the field of a single gluon of either type has been 
shown to be u T /V and uj l /V, respectively. 
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